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Introduction 

f^ The study of how photosynthetic organisms convert light offers insight not only into nature's evolutionary process, 

fS| but may also give clues as to how best to design and manipulate artificial photosynthetic systems - and also how far 

,—1 we can drive natural photosynthetic systems beyond normal operating conditions, so that they can harvest energy 

3 for us under otherwise extreme conditions. In addition to its interest from a basic scientific perspective, therefore, 
the goal to develop a deep quantitative understanding of photosynthesis offers the potential payoff of enhancing our 
' current arsenal of alternative energy sources for the future. 

i_i In the following Chapter, we consider the excitation dynamics of photosynthetic membranes in Rps. Photometricum 

cj^ purple bacteria. The first studies on purple bacteria photosynthetic membranes were concerned with the complex 

O underlying detailed structure [121 [571 1311 HS] . The interested reader might find helpful the first section of this chapter 

. where we present a summary of structures and processes that allow photosynthesis in purple bacteria. As improved 

"trt resolution became available for light-harvesting structures, so too the community's interest increased in understanding 

r^ the details of the rapid femto- to picosecond timescales for the excitation transfer process within a given harvesting 

I complex - even to the regime where quantum effects are expected. Indeed such quantum effects have recently been 

'^ confirmed in, for instance, the Fenna-Matthews-Olson complex of green sulfur bacteria [3]. However, the processes 

5 occurring at this level of detail in terms of both structure and properties of the excitation harvesting, have not yet been 

shown as being crucial to the performance of the full harvesting membrane, nor with the primary goal of any given 

photosynthetic organism: to fuel its metabolism. We focus on the transfer among different complexes, in particular the 

inter-complex excitation transfer. We consider first a system of few complexes in order to understand the consequences 

on the dynamics of the complexes' connectedness. We also consider the relative amounts of harvesting complexes 

(i.e. stoichiometry) on small sized networks in order to help establish our understanding of the behavior of complete 

^^ chromatophore vesicles. 

^-H As a whole, chromatophore vesicles comprise sections of the purple bacteria cytoplasmic membrane where a large 

^— ^ number of harvesting complexes accommodate. The current capabilities to dissect these vesicles on the nanoscale 

f — ■ through Atomic Force Microscopy (AFM) has provided evidence of changes in the conformation of chromatophores 

^^ as a result of different environmental conditions. According to Ref . [34] , membranes grown under Low Light Intensity 

^~~^ (LLI, 10 Watt/m^) present a relative amount of harvesting complexes which is different to the stoichiometry observed 

. . for bacteria grown under High Light Intensity (HLI, 100 Watt/m^), thereby pinpointing the importance of global 

^ changes in the complete vesicles as an important means of fulfilling the bacteria's metabolic requirements. In order to 

k>( study these global conformational changes - given the fact that the inter-complex transfer time-scale involves several 

; I picoseconds and the excitation delocalization length is not expected to be beyond a single harvesting complex - we 

^ choose a model of excitation dynamics which is based on a classical random walk. This random walk is coupled to 

the main processes leading to electron/chemical energy transformation, and hence the bacteria's metabolic demands. 

This Chapter is organized as follows. Section 1 provides a review of the basic structures involved in excitation 

transfer, along with a summary of the processes required for electron/chemical energy conversion. A discussion of the 

dynamics of excitations in a few model architectures is explored in section 2, in order to understand the results from 

complete LLI and HLI adapted chromatophores in section 3 and to guide our development of an analytical model in 

section 4 for determining both the efficiency and power output of any given chromatophore vesicle under arbitrary 

light intensity regimes. Lastly, in section 5 we explore the effect of incident light with extreme photon arrival statistics, 

on the resulting electronic/chemical energy conversion, in order to heuristically provide a survivability margin beyond 

which terrestrial bacteria could not survive. This is motivated by the fact that one day, it may be necessary to send 

simple bacteria into deep space and/or cope with extreme photon conditions here on Earth as a result of a catastrophic 

solar change. 
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FIG. 1: Representation of a lake containing both aerobic and anaerobic phototrophic organisms. Note that purple bacterial 
photosynthesis is restricted to the lower anaerobic layer and so they only receive solar energy that has been filtered, mainly by 
chlorophylls belonging to algae, cyanobacteria and plants. 

1. Important processes for solar /chemical energy conversion 

Purple bacteria sustain their metabolism using photosynthesis in anaerobic conditions and under the dim light 
excitation proper of several meters deeps at ponds, lagoons and streams [30j . As depicted in Figfll aerobic organisms 
are present near the surface of water reservoirs, collect blue and red spectral components of sun's light, leaving only 
the green and far red (> 750nm) components, from where purple bacteria must fulfill their energy requirements. 

The light energy absorption is accomplished through intracytoplasmic membranes where different pigment-protein 
complexes accommodate. Light Harvesting complexes (LHs) have the function of absorbing light and transfer it to 
Reaction Centers (RCs), where a charge separation process is initiated [Tj- The unpaired charge reduces a quinone, 
which using a periplasmic hydrogen, converts to quinol (Q_bH2). RCs neutrality is restablished thanks to cytochrome 
cyt charge carrier, which after undocking from the RC, must find the bci complex to receive an electron and start 
its cycle all over again. The electron in bci is given due to cytoplasmic Quinol delivery in bci. The proton gradient 
induced by the charge carriers cycling becomes the precursor of adenosyn triphosphate (ATP) synthesis: ADP+P-t> 
ATP+Energy, where ADP and P refer respectively, to adenosine diphosphate and phosphorous. The cycle is depicted 
in Figll 

1.1 Harvesting function 

Light absorption occurs through organic molecules, known as chromophores, inserted in protein complexes. Bacte- 
riochlorophyls (BChl) and carotenoids (Car) chromophores are the main absorbers in purple bacteria photosynthesis, 
principally in the far red and green, respectively. The light absorption process occurs through chromophore's Qy 
electronic transition excitation. Several chromophores are embedded in protein helices, named a and /3 apoproteins, 
inside complexes, classified by their absorption spectral maximum. Light Harvesting complex 2 (LH2) reveals two 
concentric subunits that according to their absorption maxima are called B800 or B850, composed of nine pairs 
of apoproteins |27| comprising an inner a-helix and an outer /3-helix both crossing the harvesting membrane from 
periplasm to cytoplasm, in an a/3 unit that serves to anchor a highly interacting B850 dimer and one B800 chro- 
mophore. Hence, the B800 is composed of nine chromophores, while B850 include eighteen BChla chromophores, 
having dipole moments parallel and nearly perpendicular to the membrane plane, respectively. Raman spectra using 
different excitation wavelength |13j and stoichiometry analysis [T], indicated that one carotenoid (Car) is present per 
each a/3 unit. Light Harvesting complex 1 (LHl) absorbs maximally at 883 nm, and contain 32 BChls, arranged in 16 
bi-chromophore a/3 units, surrounding an RC [22j in the same geometrical arrangement as B850 chromophores. An 
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FIG. 2: Schematic representation of the photosynthetic apparatus in the intracytoplasmic membrane of purple bacteria. The 
reaction center (RC, red) is surrounded by the light-harvesting complex 1 (LHl, green) to form the LHl+RC complex, which is 
surrounded by multiple light-harvesting complexes (LH2s) (blue), forming altogether the photosynthetic unit (PSU). Photons 
are absorbed by LHs and excitation is transferred to the RC initiating a charge (electron-hole) separation. Electrons are 
shuttled back by the cytochrome C2 charge carrier (blue) from the ubiquinone-cytochrome bci complex (yellow) to the RC. The 
electron transfer across the membrane produces a proton gradient that drives the synthesis of ATP from ADP by the ATPase 
(orange). Electron e~ is represented in blue, and quinones Qb, likely confined to the intramembrane space, in black. 



RC presents a highly interacting dimer, the special pair (P), that is ionized due to the electronic excitation transferred 
from the surrounding LHl. 

The complexes' photon absorption cross section has been calculated for LHl and LH2 complexes, where all absorbing 
molecules and extinction coefficients [T2] have been taken into account. A photon of wavelentht A, is part of the power 
spectrum of a source with occupation numbers n(A). Normalized to 18 W/m^ intensity, the rate of photon absorption 
for circular LHl complexes in Rb. sphaeroides |14j : 



7f 



n{X)a-^^l{X)dX ^ 18s- 



(1) 



The same procedure applied to LH2 complexes, yields a photon capture rate of 7^ ~ 10s~^. Since these rates are 
normalized to 18 W/m^, the extension to arbitrary light intensity / is straightforward. The rate of photon absorption 
normalized to 1 W/m^ intensity, will be 71(2) — l(0.55)s^^. From now on, subindexes 1 and 2 relate to quantities 
of LHl and LH2 complexes, respectively. Vesicles containing several hundreds of complexes will have an absorption 
rate: 



7A = /(7i^i + 72^2) 
where A^i(2) is the number of LHl (LH2) complexes in the vesicle. 



(2) 



1.2 Excitation transfer 



Excitation transfer happens through Coulomb interaction of electrons, excited to the Qy electronic transition in 
chromophores. The interaction energy can be formally written |41| : 



^«J = 2 ^ ^{'l>m<f>7l\V\<j>p(l)q)c+^^C-^^,CqaCno 

m,7i,p,q a,a' 



(3) 



where c^o-i '^na' are fermion creation and annihilation operators of electrons with spin a and a' , in the mutually 
orthogonal atomic orbitals (/i,„ and (/)„. The overlap {(f>m4'n\V\(f>p(f>q) is the Coulomb integral: 

{(t)m(l)p\V\(l3n(t)q) ^ dfi (1^2 CiC^l " 'r'i)<Pp{ri ~ r^)— — — (/);(rl - ■r'j)(t)q{r2 ~ r}) (4) 

that accounts on inter-molecular exchange contribution when donor and acceptor are at a distance comparable to 
the extent of the molecules, and the direct Coulomb contribution for an electron that makes a transition between 
(f)rn and (/)p, having both a finite value near the position of the donnor rl ~ rp, while another electron is excited 
between (/>„ and (j)q at the acceptor coordinate i^ ~ ta- In this latter situation, a commonly used framework concerns 
a tight-binding Hamiltonian, where details in the specific molecular orbitals |0„) involving mainly the Qy orbitals is 
left aside, and emphasis relies on occupation with a single index labeling an electronic state \i) concerning occupation 
in a given chromophore. The Hamiltonian H in the chromophore site basis |i), 

H = eY,\^){^\+Y.V^,\^)(J\ (5) 

i i,j 

has diagonal elements {i\H\i) — e, concerning the energy of the excitation, usually measured from the ground electronic 
state. Neighboring chromophores are too close to neglect their charge distribution and its interaction is determined 
such that the effective Hamiltonian spectrum matches the spectrum of an extensive quantum chemistry calculation 
|18j . In the LH2 complex, the B850 ring, with nearest neighbors coupling {i\H\i + 1) = 806 cm~^ or 377 cm ^^ for 
chromophores respectively, within or in different neighboring a/3 units. For next-to-neighboring chromophores the 
dipole-dipole approximation is usually used, 

with i ^ j, i = j j^ ±1 (6) 

where fli is the dipole moment and r is the distance between the interacting dipoles with constants: e = 13059 cm~^ 
and C = 519 044 A^ cm~^. The LHl complex has a spectrum maximum near 875 nm, with inter-complex distances and 
nearest neighbor interactions equal to the ones provided in the B850 LH2 ring. The molecular nature of chromophores 
involve vibrational degrees of freedom that provide a manifold within each electronic states which should be accounted. 
However, given these parameters within a harvesting complex, intra-complex energy transfer involves sub-picosecond 
time-scales, that imply relevance of electronic quantum features over the influence of thermalization in the vibrational 
manifold. If decoherence sources are not important within the time-scale of excitation dynamics on donor or acceptor 
aggregates, they can be known from the eigenvector problem H\(j)k) = Ek\4>k) whose solution provides excitonic fully 
delocalized states 

\(t>k) = ^akAi) (7) 

n 

considered to best describe the B850 ring. A slightly less interaction strength allows thermalization through the 
vibrational manifold in comparable time-scales, degrading delocalization in the B800 ring over 2-3 pigments, however, 
able to improve robustness of B800— J-B850 energy transfer |2D! . 

The fact that states involved in excitation transfer in molecules include the vibrational manifold, makes enormous the 
Hilbert space over which the sums of eq. ([3]) should be formally performed. When the transfer occurs with chromophores 
that belong to different complexes, say inter-complex energy transfer, rapid thermalization out stands direct Coulomb 
mechanism, and no defined phase relationship between donor and acceptor electronic states is expected. The hermitian 
nature of excitonic exchange is replaced due to decoherence, by a rate to describe electronic excitation transfer within 
thermalization in the vibrational manifold. According to Fermi's golden rule, adapted to the vibrational continuum, 
the rate of transfer is given by [17] : 

koA -'-^[dEf dwA I rf^^ g^("l>) expj-^VfeBT) gAJn^A) exp(-^^/fcBT) ^^^^^, 

(8) 

interpreted as a the sum of Coulomb contributions from electrons at donnor (D) and acceptor (A) aggregates in 
ground or excited states (these latter dennoted by *), Uda = (^_D*^yi|^DA|*D^A*), weighted by Boltzmann factor 
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FIG. 3: Energy level scheme of donor and acceptor molecules. Although the zero phonon line might be difTerent between both, 
energy conservation on the transfer applies due to vibrational levels. 



(kg is Botlzmann constant) and vibrational manifold multiplicity igD(A)iwD(A)) for donnor (acceptor)) at the elec- 
tronic energies Wjjim measured from the donnor (acceptor) zero phonon lines E£,ij^\^ . With the Born-Oppenheimer 
approximation [T71I1T] |^) is assumed as products of electronic |0) and vibrational \x) molecular states: 



UdA ~ {(j)D-'(l)A\VDA\(l)D(liA-) y^ {x{w*d)\x{wd)){x{wa)\x{w*a)) 
~ UDA{xi.W*D)\x{wD)){x{wA)\x{w*A)) 



(9) 



where Uda = ('/>d*'0/i|^da|''/'dV'a*)- Using the above approximation, the expression ([8]) now including the overlap 
between vibrational levels can be cast in a more illustrative form: 
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(10) 



Here, Gd{E) and Ga{E) are often called the Franck-Condon weighted and thermally averaged combined density of 
states. Explicitly: 



Gd{E) 



dw'\ 



gl,{wl,)ciqi{-w*^/kBT)\{x{w*o)\x{wD))\" 



with equal expression for acceptor molecule by replacement D —^ A. Forster showed [TT] that these distributions are 
related to extinction coefficient e{E) and fluorescence spectrum foiE) of direct experimental verification: 
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where Nq — 6.022 x 10^*^ is the number of molecules per mol per cm'^, n is the refractive index of the molecule sample, 
c the speed of light, and tq the mean fluorescence time of the donnor excited state. For normalized spectra 
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(12) 



and from the relations wu = E'jj^, + w'^ — E and w'^ — —Eaq + wa + E schematically presented in FigjSl the Forster 
rate is cast 



... = |ia„.i^/.^H<S^. 



(13) 



Therefore, whenever fluorescence and absorption spectra are available, an estimate for the excitation transfer rate can 
be calculated. 

Thermalization occurs flrstly in the vibrational manifold of the electronic states involved. Due to the greater energy 
gap of electronic transitions compared with the one of vibrational nature, on a longer time-scale thermalization also 
occurs in excitonic states. Accordingly, the calculation in eq.(|9| involves an statistical thermal mixture, explicitly: 



P = 



^^exp(-Sfc//cBT)|,/)fe)(0,| 



Tr{.} 



(14) 



where Tr{-} is trace of the numerator operator, used to normahze the state. Hence, in a straightforward fashion 

1 



Uda = Tt:{pVda} = 



1 

j:^exp{-Ek/kBT) 



J2 J2 cM-Ek/kBT){cl)k\VDA\(bp) 



k V 



/ . ak,iap^jVij 



(15) 
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where the element {(j)k\yDA\<t>p) are the elements of interaction among exciton states in molecules on different complexes 
and use is made of the individual contributions of excitonic states in chromophore's site basis eq.([7|. 

Summarizing, excitation transfer occurs through induced dipole transfer, among BChls transitions. The common 
inter-complex BChl distances 20-100 A [2j [34] cause excitation transfer to arise through the Coulomb interaction on 
the picosecond time-scale }17j . while vibrational dephasing destroys coherences within a few hundred femtoseconds 
[531 US]- As noted, the Coulomb interaction as dephasing occurs, makes the donor and acceptor phase become uncor- 
related pointing into a classical rate behavior. Transfer rate measures from pump-probe experiments agree with the 
just outlined generalized Forster calculated rates |17) , assuming intra-complex delocalization along thermodynamical 
equilibrium. LH2— >-LH2 transfer has not been measured experimentally, although an estimate of ^22 = 10 ps has been 
calculated [17]. LH2— >■ LHl transfer has been measured for R. Sphaeroides as ^21 = 3.3ps 16 . Due to formation 
of excitonic states [21], back-transfer LHl— >■ LH2 is enhanced as compared to the canonical equilibrium rate for a 
two-level system, up to a value of ti2 = 15.5ps. The LHl— >-LHl mean transfer time in has not been measured, but 
the just mentioned generalized Forster calculation [3T] has reported an estimated mean time in of 20 ps. LHl— >■ 
RC transfer occurs due to ring symmetry breaking through optically forbidden (within ring symmetry) second and 
third lowest exciton lying states [18], as suggested by agreement with the experimental transfer time of 35-37 ps at 
77 K [43l|44]. Increased spectral overlap at room temperature improves the transfer time to ti^Rc = 25 ps [42]. A 
photo-protective design makes the back-transfer from an RC's fully populated lowest exciton state to higher-lying 
LHl states occur in a calculated time of taci =8.1 ps [TH], close to the experimentally measured 7-9 ps estimated 
from decay kinetics after RC excitation [ID]. 

Table IT] shows the results of mean transfer times presented in Ref. [31] through the above mentioned calculation, 
compared with the experimental evidence restricted to different complex kind from the spectral resolution requirement 
of pump-probe spectroscopy. Since LHl<->- LHl and LH2o LH2 transfer steps involve equal energy transitions, no 
experimental evidence is available regarding the rate at which these transitions occur. The experimentally determined 
B800— > B850 rate was l/700fs 38J. The inter-complex transfer rate between LH2-H- LHl have been determined 
experimentally to be 1/3. 3ps ^]. Experimentally, LH1<-^ RC forward transfer rate ranges between l/50ps and 
l/35ps, while back-transfer rate ranges between l/12ps and l/8ps [39l|40l[43]. It is interesting to note that exists a 
two fold difference in the experimental and theoretical determined LH2— >-LHl, ascribed to BChla Qy dipole moment 
underestimation. It is assumed for theoretical calculation a value of 6.3 Debye, while a greater BChla Qy dipole 
moment in Prostecochloris aestuarii (not a purple bacterium) of 7.7 Debye has been determined [311 . On the other 
hand, LHl— >RC theoretical calculation gives a greater value for tranfer rate, thought to arise due to an overestimate 
of LHl exciton delocalization [31]. This rate decreases when delocalization is assumed over fewer BChl's, therefore, 
further research is needed to understand the effect of decoherence sources (static inhomogeneities and dynamical 
disorder due to thermal fluctuations) on the delocalization length. 
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TABLE I: Theoretical estimation/experimental evidence of inter-complex transfer times in picoseconds. N.A are not available 
data. 



1.3 Dissipation 



Excitation in chromophores might be dissipated by two main mechanisms. The first is fluorescence, where the elec- 
tronic excited state has a finite lifetime on the nanosecond time-scale, due to its interaction with the electromagnetic 
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FIG. 4: Dissipation mechanisms. In (a), the electron de-excites due to its interaction with the quantized electromagnetic 
vacuum field through a fluorescent photon. In (b) internal conversion mechanism, where the vibrational levels overlap induces 
a transition between electronic excited and ground states. Dissipation overcomes when thermal equilibrium is reached in the 
vibrational manifold of electronic ground state. 

vacuum |37| . The second is internal conversion, where the electronic energy is transferred to vibrational degrees of 
freedom. 

Within the Born-Oppenheimer approximation, the molecular state ^, can be decomposed into purely electronic 
(j> and (nuclear) vibrational x states. The transition probability between initial state ^^ and final state ^I'/, is 
proportional to {'i'i\H\'$f) ex (xilXf)- Note that Xi iXf) ^-^e vibrational levels in the ground (excited) electronic state 
manifold (see FigHl). If the energy difference is small, and the overlap between vibrational levels of different electronic 
states is appreciable, the excitation can be transferred from the excited electronic state, to an excited vibrational level 
in the ground electronic state. This overlap increases with decreasing energy difference between electronic states. As 
higher electronic levels have smaller energy difference among their zero phonon lines, internal conversion process is 
more probable the higher energy electronic states have. Fluorescence and internal conversion between first excited 
singlet and ground electronic states, induce dissipation in a range of hundreds of picoseconds and a few nanoseconds. 
Numerical simulations are performed with a dissipation time including both fluorescence and internal conversion of 
I/td = 1 ns, also used in J3T]. 

1.4 Special pair (SP) ionization 

From the LHl complex excitation reaches the RC, specifically the special pair (SP) dimer. The excitation can be 
transferred back to its surrounding LHl, or initiate a chain of ionizations along the A branch, probably, due to a 
tyrosine residue strategically positioned instead of a phenylalanine present in the B branch [531. Once the special 
pair is excited, it has been determined experimentally lOj that takes 3-4 ps for the special pair to ionize and produce 
a reduced bacteriopheophytin, H^, in a reaction SP* — > SP^H^. This reaction initiates an electron hop, to a 
quinone Qa in about 200 ps, and to a second quinone, Qb if available. Initially, the ionized quinol Qj^ captures an 
introcytoplasmic proton and produces hydroxiquinol QbH, which after a second ionization that produces QbH^ to 
form quinol Q_b^2- After any SP ionization a neutrality restablishment is required, provided by the cytochrome cyt 
charge carrier. After SP ionization, the cytochrome diffuses from the bcl complex to a RC in order to replenish its 
neutrality SP"*" — )-SP, within several microseconds [H]. The first electron transfer step P* — >■ P"*" occurs in the RC 
within t+ =3 ps, used for quinol {QbH2) production |17j . 

1.5 Quinone-quinol cycling 

The RC cycling dynamics also involves undocking of QbH2 from the RC due to lower affinity among RC and this 
new product. Quinol starts a migration to the hcl complex where enables the ionization of the cytochrome cyt charge 
carrier, while a new quinone Qb molecule docks into the RC. The time before quinol unbinds, and a new Qb is 
available, has been reported within milliseconds i28 to highlight quinol removal o as the rate limiting step [3H] if 
compared to special pair restablishment. 

Even though it has been reported that excitation dynamics change as a function of the RCs state ^ 18j , at a first 
glance the several orders of magnitude difference among the picosecond transfer, the nanosecond dissipation and the 




FIG. 5: Schematic of the biomolecular photosynthetic machinery in purple bacteria, together with relevant inter-complex mean 
transfer times tij, dissipation rate 7d, and normalized light intensity rate 71(2) 

millisecond RC cycling, seems to disregard important effects due to these mechanisms' interplay. However, the quinol- 
quinone dynamics leaves the RC unable to promote further quinol production and eventually enhances the influence 
of dissipation of a wandering excitation, evident when none RC is available and the unique fate of any excitation is 
to be dissipated. 

Interestingly, the quinone- quinol mechanism has been well established and thought to be of priority on adaptations 
of bacteria, that seem to respond to its dynamics. For instance, an observed trend for membranes to form clusters 
of same complex type [33j seems to affect diffusion of quinones, enhanced when, due to higher mobility of LHls, 
left void spaces help travel quinones to the periplasm. Negligible mobility of LH2s in their domains, would restrict 
metabolically active quinones to LHl domains [35|. Easier diffusion of quinones, quinol and cytochromes promotes 
higher availability of charge carriers in RC domains under LLI conditions, increasing the rate at which RCs can 
cycle. The RC cycling dynamics and its connection to the membranes performance has been accounted in [51 [S] in a 
quantitative calculation to understand the effect of core clustering and stoichiometry variation in the RC supply or 
in the efficiency of the membranes from experimentally obtained Atomic Force Microscopy images, to be presented in 
this chapter. 

2. Exciton kinetics 

Figure p] summarizes the relevant biomolecular complexes in purple bacteria Rsp. Photometricum [33 , together 
with experimental- theoretical if the former are not available- timescales governing the excitation kinetics: absorption 
and transfer; and reaction center dynamics: quinol removal. 



2.1 Model 



The theoretical framework used to describe the excitation transfer must be built around the experimental (if 
available) and theoretical parameters just outlined. Remind that the thermalization process occurs faster than inter- 
complex energy transfer, and provides the support to rely in a classical hopping process, since phase information is 
lost well within the time frame implied by direct Coulomb coupling. Accordingly, we base our analysis on a classical 
random walk for excitation dynamics along the full vesicle, by considering a collective state with N = A'2 + 2A^i sites - 



9 

resulting from N2 LH2s, A^i LHls and hence Ni RC complexes in the vesicle. The state vector p= {pi,p2, ■■■,Pm) has 
in each element the probability of occupation of a collective state comprising several excitations. If a single excitation 
is allowed in each complex, both excited and ground states of any complex should be accounted and the state space 
size is Af = 2 X 2 X 2... = 2^. On the other hand, if only one excitation that wanders in the whole network of 

N 

complexes is allowed, a site basis can be used where each element of the state vector gives the probability of residence 
in the respective complex, and reduces the state vector size to M — N . In either case the state vector time evolution 
obeys a master equation 

M 

dtpm = Y.^^,,pAt)- (16) 

where Gi^ is the transition rate from a collective state or site i- whether many or a single excitation are accounted, 
respectively - to another collective state or site j. Since the transfer rates do not depend on time, this yields a 
formal solution p{t) = e'^*p{0). However, the required framework depends on exciton abundance within the whole 
chromatophore at the regime of interest. 

For instance, purple bacteria ecosystem concerns several meters depths, and should be reminded as a low light 
intensity environment. Within a typical range of 10-100 W/m^ and a commonly sized chromatophore having sa 400 
LH complexes, eq.([2| leads to an absorption rate ja ~100-1000 s~^, which compared with the dissipation mechanisms 
(rates of w 10^ s^^) imply that an absorption event occurs and then the excitation will be trapped by a RC or become 
dissipated within a nanosecond, and other excitation will visit the membrane not before some milliseconds have 
elapsed. However, it is important to remind the nature of thermal light where the possibility of having bunched small 
or long inter-photon times is greater than evenly spread, with greater deviations from poissonian statistics the grater 
its mean intensity is. Therefore, regardless of such deviations, under the biological light intensity conditions, the 
event of two excitations present simultaneously along the membrane will rarely occur and a single excitation model 
is accurate. 

2.2 Small architectures 

Small absorption rates lead to single excitation dynamics in the whole membrane, reducing the size of p{t) to the 
total number of sites N. The probability to have one excitation at a given complex initially, is proportional to its 
absorption cross section, and can be written as p{Q) — — (71, ...,72, ..., 0,..), where subsets correspond to the A''i 

Ni N2 Ni 

LHls, the N2 LH2s and the A^i RCs respectively. 

2.2.1 Complexes arrangement: architecture 

To gain physical insight on the global behavior of the harvesting membrane, our interest lies in the probability to 
have an excitation at a given complex kind k € LH1,LH2 or RC, namely pk, given that at least one excitation resides 
in the network: 

- f-^\ Pf^^^^ (^'7\ 

Pkit) = —^ — . (17) 

The effects that network architecture might have on the model's dynamics, are studied with different arrangements 
of complexes in small model networks, focusing on architectures which have the same amount of LHl, LH2 and RCs 
as shown in the top panel of FiglHFa), (b) and (c). The bottom panel Fig|6] (d)-(e)-(f) shows that pk values for RC, 
LHl and LH2 complexes, respectively. First, it is important to notice that excitations trend is to stay within LHl 
complexes, and not in the RC. Fig|6][d) shows that the highest RC population is obtained in configuration (c), followed 
by configuration (a) and (b) whose ordering relies in the connectedness of LHls to antenna complexes. Clustering of 
LHls will limit the number of links to LH2 complexes, and reduce the probability of RC ionization. For completeness, 
the probability of occupation in LHl and LH2 complexes (FigslG^e) and (f), respectively), shows that increased RC 
occupation benefits from population imbalance between LHl enhancement and LH2 reduction. As connections among 
antenna complexes become more favored, the probability of finding an excitation on antenna complexes will become 
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FIG. 6: Top panel: Three example small network architectures. The bottom panel shows the normalized probabilities for finding 
an excitation at an RC (see (d)), an LHl (see (e)), or an LH2 (see (f)). In panels (d)-(f), we represent these architectures as 
follows: (a) is a continuous line; (b) is a dotted line; (c) is a dashed line. 

smaller, while the probability of finding excitations in RCs is enhanced. This prchminary result, illustrates that if 
the apparent requirement to funnel excitations down to RCs in bacterium were of primary importance, the greatest 
connectedness of LH1-LH2 complexes should occur in nature as a consequence of millions of years evolution. However, 
as will be presented, the real trend to form LHl clusters, reduces its connectedness to antenna LH2 complexes and 
somehow pinpoints other mechanisms as the rulers of harvesting membranes conformation and architecture. 



2.2.2 Relative amount of complexes: Stoichiom,etry 

We can also address with use of small architectures the effect of variation in the relative amount of LH1/LH2 
complexes, able to change the population of the available states. FigjT] shows small networks of LH-RC nodes, where 
the relative amount of LH2 and LHl complexes quantified by stoichiometry s — N2/N1 is varied, in order to study 
the exciton dynamics. In FigJ8[a) the population ratio at stationary state of LHs demonstrate that as stoichiometry 
s becomes greater, the population of LHls, becomes smaller, since their amount is reduced. It is apparent that RC 
population is quite small, and although their abundance increases the exciton trend to be found in any RC (FiglSTb)), 
generally, excitations will be found in harvesting complexes. The population of LHs should be dependent on the ratio 
of complexes type. As verified in FigMb), RCs have almost no population, and for the discussion below, they will 
not be taken into account. Populations can be written as: 



hit 

P2it 



00) = /i(s) 

00) = /2(s) 



iVi 



Ni +N2 
N2 

N1+N2 



l + s 

l+S 



where the dependence on the amount of complexes is made explicit with the ratio 



JVfc_ 



(18) 
(19) 
and where /i(s) and /2(s) 



are enhancement factors. This factor provides information on how the population on individual complexes changes. 
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FIG. 7: Networks with diflterent stoichiometries, from left to right, top to bottom, s={1.04, 2.06, 3.08, 4.44, 5.125, 6, 7.16, 
11.25, 15.33, 23.5, 48}, and equal number of harvesting complexes. 



beyond the features arising from their relative abundance. With use of eqs.(18p9), /i(s) and /2(s) can be numerically 



calculated provided that pkit — >■ cxd) can be known from the master equation, while s is a parameter given for each 
network. The results for enhancement factors are presented in FigMc). The enhancement factor /2(s) for LH2 seems 
to saturate at values below one, as a consequence of the trend of excitations to remain in LHls. This means that 
increasing further the number of LH2s will not enhance further the individual LH2 populations. On the other hand 
/i(s) has a broader range, and increases with s. This result reflects the fact that population of individual LHls will 
become greater as more LH2 complexes surround a given LHl. An unconventional architecture (third column, second 
row in FigJTl) has an outermost line of LHl complexes, whose connectedness to LH2s is compromised. In all the results 
in Fig. Is] (sixth point), this architecture does not follow the trends just pointed out, as LHl and RC population, 
and enhancement factors, are clearly reduced. The population of LHl complexes depends on their neighborhood 
and connectedness. Whenever connectedness of LHl complexes is lowered, their population will also be reduced. 
Hence, deviations from populations trend with variation of stoichiometry, are a consequence of different degrees of 
connectedness of LHls. 

Up to this point, the master equation approach has helped us understand generally the effect of stoichiometry and 
architecture in small networks. Two conclusions can be made: 

1. Connectedness of LH2 complexes to LHls, facilitates transfer to RCs 

2. The relative amount of LH2/LH1 complexes, namely, stoichiometry s = N2/N1, when augmented, induces 
smaller population on LHl-RC complexes. On the other hand, smaller s tends to increase the connectedness of 
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FIG. 8: In (a) stationary state populations for LII2s (circles), LHls (diamonds) and RCs (crosses), as a function of the 
stoichiometry of membranes presented in FiglT] In (b) a zoom of RC populations is made, and in (c) the enhancement factors 
/i(s) (diamonds) and /2(s) (circles) are presented. 
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FIG. 9: Normalized probabilities p^ for finding the excitation at an LII2 (dashed), LHl (dotted) or at an RC (continuous), for 
(a) t+ — 3ps, and (b) i+ — !> cxj. Crosses are the results from the Monte Carlo simulation. 



LHls to LH2s and hence, the population of individual LHl complexes. 



2.2.3 Special pair ionization 



Another basic process involved in the solar energy conversion is the ionization of the special pair in the RC, and 
eventual quinol QbH2 formation. Remind that once quinol is formed, the special pair is unable to use further incoming 
excitations before quinol undocks and a new quinone replaces it. Even though the RC neutrality-diifusion process 
is propelled by complicated dynamics and involved mechanisms, in an easy approach, let us assume that the RC 
dynamics will proceed through a dichotomic process of "open" and "closed" RC states. In the open state, special pair 
oxidation is possible, while when closed, special pair oxidation to form quinol never happens, hence t-|_ — >• oo 

The eifect of open and closed RC states changes the exciton kinetics. We start with a minimal configuration 
corresponding to a basic photosynthetic unit: one LH2, one LHl and its RC. FigureloFa) shows that if the RC is open, 
excitations will mostly be found in the LHl complex, followed by occurrences at the LH2 and lastly at the RC. On 
the other hand. Figure lO^b)) shows clearly the different excitation kinetics which arise when the RC is initially unable 
to start the electron transfer P* -^ P"*", and then after « 15ps the RC population becomes greater with respect to the 
LH2's. This confirms that a faithful description of the actual photosynhesis mechanism, even at the level of the minimal 
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unit, must resort into RC cycling, given that its effects are by no means negligible. Moreover, comparison among 
FigsJ6[d) andB^also presents a feature that is usually undermined when small architectures are used to straightforward 
interpret its results as truth for greater, real biological vesicles. Energy funneling becomes smaller with the number 
of antenna LH2 complexes, thereby, in architectures with many harvesting antenna complexes, excitation will find 
it more difficult to arrive to any of the relatively spread RCs. Besides, although LH2— >-LHl transfer rate is five-fold 
the back-transfer rate, the amount of smaller sized LH2s neighboring a given LHl will increase the net back-transfer 
rate due to site availability. Hence, the funneling concept might be valid for small networks [171 131] . however, in 
natural scenarios involving entire chromatophores with many complexes, energy funneling might not be priority due 
to increased number of available states, provided from all LH2s surrounding a core complex, and globally, from the 
relative low RC abundance within a real vesicle. 

It is important to mention that results for master equation calculations require several minutes in a standard 
computer to yield the results shown in Fig|8J and that these networks have an amount of nodes an order of magnitude 
smaller than the actual chroniatophore vesicles. Dynamics concerning the RC cycling have not been described yet, 
fact that would increase further the dimension of possible membrane's states. To circumvent this problem, further 
analysis will proceed from stochastic simulations^ and observables will be obtained from ensemble averages. 

2.3 Full vesicles 

A real vesicle involves several hundreds of harvesting complexes. Given the large state-space needed to describe 
such amount of complexes and our interest to inquire on a variety of incoming light statistics in the sections ahead, our 
subsequent model analysis will be based on a discrete-time random walk for excitation hopping between neighboring 
complexes. 

Simulation algorithm. 

In particular, we use a Monte Carlo method to simulate the events of excitation transfer, the photon absorption, the 
dissipation, and the RC electron transfer. We have checked that our Monte Carlo simulations accurately reproduce 
the results of the population-based calculations described above, as can be seen from FigsMa) and (b). The Monte 
Carlo simulations proceed as follows. In general, any distribution of light might be used with the restriction of having 
a mean inter-photon time of 7^ from eq. i2h . Accordingly, a first photon is captured by the membrane and the time 
for the next absorption is set by inverting the cumulative distribution function from a [0,1] uniformly distributed 
(Unit Uniformly Distributed, UUD) random number. This inversing procedure is used for any transfer, dissipation 
or quinol removal event as well. The chosen absorbing complex is randomly selected first among LHl or LH2 by 
a second UUD number compared to the probability of absorption in such complex kind, say A'^x(2)7i(2)/7A i^i LHl 
(LH2), and a third UUD random number to specifically select any of the given complexes, with probability 1/Ni(^2)- 
Once the excitation is within a given complex, the conditional master equation given that full knowledge of the 
excitation residing in site z, only involves transfers outside such site, say dtPi = — (X^i l/^i,i + lD)Pi, whose solution 
is straightforward to provide the survival probability and its inverse, of use to choose the time t* for the next event 



according to eq.(16l from a UUD number r: — logr/(^ ■ l/t^ j + 7d) = i* ■ Once t* is found, a particular event 
is chosen: transfer to a given neighboring complex j with probability (l/ii.j)/(X)7 ^/^i-i ^ ^d) or dissipation with 
probability Jd/C^j ^/ti,j +7£i), which are assigned a proportional segment within [0,1] and compared with another 
UUD number to pinpoint the particular event. If the chosen event is a transfer step, then the excitation jumps to the 
chosen complex and the transfer-dissipation algorithm starts again. If dissipation occurs, the absorption algorithm is 
called to initiate a new excitation history. In a RC, the channel of quinol ionization is present with a rate l/t+ in 
an event that if chosen, produces the same effect as dissipation. Nonetheless, the number of excitations that become 
SP ionizations are counted on each RC, such that when two excitations ionize a given RC and produce quinol, it 
becomes closed by temporally setting l/t+ = at such RC. Quinol unbinding will set "open" the RC, not before 
the RC-cycling time with mean t, has elapsed, chosen according to a poissonian distribution. The algorithm can be 
summarized as follows: 

1. Create the network: Obtain coordinates and type of LHs, and label complexes, for instance, by solely numerating 
them along its type, say complex 132 is of type 2 (we use 1 for LHl, 2 for LH2 and 3 for RC). Choose the j 
neighbors of complex i according to a maximum center to center distance less than ri + r2 + S , r2 + r2 + S and 
ri+ri + S for respective complexes. We use S — 20 A, chosen such that only nearest neighbors are accounted and 
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further increase of S makes no difference on the amount of nearest neighbor connections, although further increase 
may include non-physical next to near-neighbors. In practice, the network creation was done by three arrays, 
one, say neigh{i,j) with size M x S, with M complexes as described above, and S as the maximum number of 
neighboring complexes among all the sites, hence requiring several attempts to be determined. Minimally j < I 
for an LH2, concerning the dissipation channel, j < 2 for LHl including both dissipation and transfer to its RC, 
and J < 3 for a RC accounting on dissipation, RC ionization an transfer to its surrounding LHl. The other 
arrays are built, say size{i), with M positions, that keep on each the number of neighbors of the respective i 
labelled complex, and rates(i,j) where at each position the inter-complex rate i — > j is saved. For instance, 
rates{i, 1) of any RC will be the ionization rate l/t+. 

2. Send photons to the network: On a time t* = —log{r)/jA according to eq.(l2|, with r being an UUD number. 
Choose an LH2 or an LHl, according to the probability of absorption from the cross section of complex type 
^i(2)li{2) /lA- Add one excitation to the network, say n = n + 1, and assign the initial position pos{n) = i of 
the excitation according to another UUD that selects an specific labelled i complex. Remind that n is bounded 
by the maximum amount of excitations allowed to be at the same time within the membrane, usually being one. 

3. If the zth complex is excited, the construction of the above mentioned arrays make the cycle of excitation 
dynamics straightforward since the network is created only once, and dynamics only require to save the complex 
i where the excitation is, and then go through cycles of size size{n, i) to acknowledge the stochastically generated 
next time for a given event. Excitation can be transferred to the available neighbors, become dissipated or a 
RC ionization event. Order all times for next events in order to know which will be the next in the array, say 
listimes{p) with p < n, where ij^jj^ — listemp{l) . In parallel, update an organized array that saves the next 
process with the number of the neighbor to which hopping occurs, or say a negative number for RC ioinization 
and another negative number for dissipation. 

4. Jump to next event: By cycling over the n present excitations, increase time up to the next event ijj^jjj. If RC 
cycling is accounted, check which time among tmin and the next opening RC time tjic (its algorithm is to be 
discussed in the following) is the closest, and jump to it. 

5. Change state of excitations or that of RCs: Update the current site of the excitation n, or whether it becomes 
a dissipation or a RC ionization. If the latter process occurs, keep in an array, say rcstate{k) whose size 
equals the total amount of RCs, the number of excitations that have become ionizations from the last time 
the fcth RC was opened. If rcstate{k) = 2 then the fcth RC is closed by redefinition of rate{i, 1) = and a 
poissonian stochastically generated opening time with mean r is generated. This time interval is kept in an 
array rctimes{k). Now, introduce this time interval into an ascending ordered list among all closed RCs opening 
times such that the minimum tuc is obtained. If tuc < ^niin then jump to that time and open the fcth RC by 
letting rcstate{k) — and ratesii, 1) = l/t+. 

6. Look which is minimum among t* , tfic and ijj^j^ and jump to steps pi or H according to whether i* < {t^Ci ^niin) 
or (iflcimin) < **' respectively. 

7. If the maximum amount of excitations chosen from the initialization, have been sent to the membrane, finish 
all processes and write external files. 

The language used to program this algorithm was FORTRAN??, to point out that these calculations do not require 
any high-level language. 

Excitation dynamics trends in many node- complexes networks. 

In order to understand at a qualitative degree the excitation dynamics trends involved in full network chro- 
matophores, a few toy architectures have been studied, shown in Fig |10| In this preliminary study it is of interest 
to understand the excitation kinetics in complete chromatophores. In particular, it is useful to understand if any 
important feature arises according to nature's found tendency of forming clusters of the same complex type. In AFM 
images |34) it has been found that there is an apparent trend to form clusters of LHl complexes with simultaneous 
formation of LH2 para-chrystalline domains. The reason that has been argued for this trend [34] involves the RC 
cycling dynamics and can be explained as follows. 

The charge carrier quinone-quinol require diffusing through the intracytoplasmic membrane within the void spaces 
left among harvesting complexes, in order to reach the bcl complex and complete the electric charge cycle with 
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cytochorme that however, diffuses through the periplasm as schematicahy shown in FigJ2J The closeness among LH2 
complexes in these para-chrystalline domains restricts the void spaces required for diffusion of quinone-quinol to the 
LHl domains, where charge separation is taking place. Then, such aggregation indeed helps to improve the time 
it takes to quinone perform the whole RC-bcl-RC cycle, by restricting its presence to RC domains. However, an 
advantage concerning excitation dynamics has been heuristically proposed, where the path an excitation has to travel 
to reach an open RC once it encounters a closed RC is reduced due to LHl clustering. In this section we investigate 
this latter possibility from a more quantitative point of view and address the former in the next section. 

In order to understand exciton kinetics, for instance, let us fix our attention on the probability of a RC ionization 
when excitations start in a given LH2 on the configurations shown in FigjTOJ To that end, let us introduce first the 
number njicij of excitations that are absorbed at site i and become ionization of the special pair at a given RC, 
labeled j. Also let f^ be the unitary vector pointing in the direction from LH2 site i to RC site j. As an analogous to 
a force field, let vector field v correspond to a weighted sum of the directions on which RC ionizations occur, starting 
for a given LH2 complex 



Y.jnBcJ^3 



E, 



n 



2 



(20) 



of help as a purely numerical calculation of graphical interpretation. The normalization gives a basis to compare 
only the directions for RC ionizations independent whether excitations are dissipated or not. Note that f^ has both 
positive and negative signs, hence the vector v^ will have greater magnitude as RC ionizations occur more frequently 
at a given direction, or less magnitude as ionizations occur in more spread out directions. 

Several features can be understood from this vector field. Figure [T0| shows the vector fields attracted to the LHl 
most abundant regions. The arrows illustrate the favored statistical direction in which excitations will become RC 
ionization, and they do not imply that an excitation will deterministically move in that direction if it were absorbed 
in such place. Firstly, in all cases, a tendency to point on clustered LHls stands from distant LH2 complexes, a fact 
that emphasizes the very rapid transfer rates that lead to excitation's rapid global membrane sampling. Secondly, it 
is evident that better excitation attractors are made with increasing the core centers cluster size, easily observed in 
the top panel. Hence, according to this result one may naively state that a funneling effect is apparent and strongly 
dependent on cluster size. However, on the middle panel, it can be seen that this funneling effect is vanishingly small 
at sites where some degree of symmetry is present concerning the distribution of neighboring RCs, even though if these 
latter display some clustering degree. Third, a trend is displayed where the spreading of the flux becomes smaller 
with a lesser amount of RCs. However it is important to be aware that in more detail and in a general fashion, one 
also sees that even the LH2s closest to a given RC, "feel" the presence of further RCs. This implies that a significant 
amount of excitations that reach LH2 complexes neighboring a given LHl, will better prefer hopping to neighboring 
LH2 complexes to eventually reach a distant RC: no "funneling", as usually understood. 

In a more realistic situation, the clustering of LHl complexes is just a trend, and in average any cluster is formed 
by a few 3-4 core complexes in HLI situations and somehow less for LLI membranes as shown in Fig |ll[ where 
empirical architectures are presented in complete accordance with experimental data, taken from Ref.[33]- Along 
these figures, the vector field calculated for the HLI and LLI empirical membranes is made, to highlight no clear trend 
of excitations to be transferred immediately to neighboring LHl complexes, and disregard a funneling LH2— >LH1 as 
the unique direction of excitation fiux when many complexes are accounted. 



In detail, from the reasons just provided, in Fig 11 a) it is apparent that in HLI membranes, individual LHl 
complexes are not efficient attractors of excitations. Here, not even the LH2s closest to a core center complex are 
notably attracted to such cores. On the other hand, a more uniform flux seems to be shared on clusters of LH2s 
in the LLI membrane (Figjlljb)), that points to locations where more core centers are displayed. Hence, it is more 
evident the directionality of the flux to the most LHl abundant locations that gives a statistical preference to ionize 
the closest RCs due their low global abundance in LLI membranes. The greater amount of LHl complexes in the HLI 
membrane induces no significant preference for excitations to hop to the closest lying RCs and provides evidence that 
the excitation process proceeds without a clear funneling effect. If a single excitation can become an ionization in a 
single RC, the flux is clear, but if there are many and spread, the flux becomes random, as would be expected from 
any random walk. In summary, these results allow a mental picture of the excitation transfer process: Excitations 
become absorbed and start wandering along the membrane in a random walk regardless of a funneling effect, up to 
the moment where they reach a RC and become ionized, in a process where statistical preference is present in LLI 
membranes due to the reduced amount of RC, and presenting more clearly the random wandering on HLI due to 
higher spread out possibilities of ionization. This fact can also be quantitatively investigated. 

Suppose that a given LHl is surrounded by a few LH2 complexes. If an excitation is absorbed into one of these 
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FIG. 10: Vector field v for different configurations, where the effect of core centers clustering is investigated (100000 runs 
starting in each site). 



neighboring LH2s, it will have a survival probability decaying with rate 1/^21 +n/t22 where n wold be 3-4 corresponding 
to the number of LH2s neighboring such absorbing LH2. Given the presented transfer times (remind ^21 = 3.3ps and 
^22 = lOps) both terms in the survival decay rate become almost equal for n = 3, and hence no preference will 
occur on the excitation to be transferred to the closest LHl. If LII2^LII1 transfer would actually happen, the same 
exercise can be done for the survival probability within the LHl complex, decaying with a rate mltyi + l/^i.flc 
(remind tyi = 15.5ps and t\^^c ~ 30ps) which assuming a single surrounding LII2 complex {m — 1) would give 
a preference for excitations to return to antenna complexes. Hence, if the availability of neighbors is accounted, 
the funneling LH2— >-LHl is lost, which along the enhanced back-transfer rate LHl— 7>LH2 as compared to LHl— >RC, 
provides the basis for understanding complete chromatophores as networks where actual performance depends upon 
the RC availability, since the event of ionization only depends on the probability to get to RC sites having no preference 
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FIG. 11: Vector field v for HLI (a) and LLI (b) membranes (100000 runs each site). 



to be visited within the network. Excitations sample great portions of the whole membrane in its hopping, hence 
become able to reach RCs far away from the absorption sites before dissipation overcomes. 

Dissipation itself also provides interesting features. The dissipation di measures the probability for excitations to 
dissipate at site i from the the total amount of excitations being absorbed in the membrane ha'- 



d, = 



UA 



(21) 



that has an straightforward relation with the global efficiency rj of the membrane, that accounts on the probability of 
any excitation to be used as an SP ionization. Given that excitations can only become RC ionizations or be dissipated, 
the sum over all complexes of dissipation probability will give the probability of any excitation to be dissipated, hence 
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1 



Y^d^ 



(22) 



Figure [12] shows the dissipation on the membranes shinned with the respective illumination rate of their growth 
(LLI 10 W/m^ and HLI 100 W/m^). Figure 12 'a) shows that the LLI membrane has highly dissipative clusters of 
LH2s, in contrast to the uniform dissipation in the HLI membrane (see Figjl^b)). 

This result addresses distinct features probably connected with the own requirements for bacteria when different 
light intensity is used during growth. Under LLI, bacterium might require use of all the available solar energy to 
promote its metabolism. In the numerical simulations (as expected from nature), the dissipation rate is set equal on 
any site, and therefore dissipation is only dependent on the probability for an excitation to be in a given complex. 
If more dissipation is found in some regions, it can there be supposed that excitations remain more time at such 
domains. Hence the tendency for excitations to dissipate implies they reside longer in LH2 complex domains and 
justify the view of LH2 clusters as excitation reservoirs. Although RC cycling has not been accounted yet, a given 
history might include an event where the excitation reaches closed clustered RCs and then jumps back to these LH2 
domains waiting before the RCs become available again. On the other hand, HLI membranes display an evenly 
distributed dissipation. For instance, at HLI not all excitations might be required and dissipation can be used as a 
photo-protective mechanism. In this case, if dissipation were highly concentrated, vibrational recombination would 
overheat such highly dissipating domains. 

Beyond these local details, an analysis regarding average values of dissipation D}. and residence pn^ probabilities 
on a complex type /c (fc = 1, 2 are respectively LHl, LH2 complexes) better supports the view of a completely random 
excitation hopping process. Table ITT] shows the numerical results of pa^ , that concerns the probability to find an 
excitation in a given complex type and is calculated from the sum of the residence times of excitations at a given 



18 




(b) 




FIG. 12: Contour plots for dissipation in LLI (a) and HLI (b) membranes. Greater contrast means higher values for dissipation. 
The simulation is shown after 10® excitations were absorbed by the membrane with rate 7a . 

complex type. It is straightforward to see that pn^ is closely related to the probability of dissipation at such complex 
type, therefore, dissipation can correctly measure where excitations are to be found in general. The randomness 
of excitation dynamics is illustrated by realizing that dissipation in a given complex type depends primarily on its 
relative abundance, since ^ « ^, and justifies that apart from LII2 clusters local variations, in the mean, all sites 
behave equivalently, and no dynamical feature arises to set a difference among LHls and LH2s able to argue on the 
stoichiometry adaptation experimentally found. Lastly, as expected from the dynamical equivalence of sites, along the 
result obtained with toy architectures of varying stoichiometry where RC population heavily rose due their abundance 
but relatively invariant to arrangement, the probability to reach a RC solely depends on its abundance, and therefore 
the global efficiency ry is greater in HLI membranes. 



Membrane 


PR2 


PRi 


Z>2 


Di 


D2 


„_ iV2 
*- iVi 


r? = l-^ 


LLI 
HLI 


0.72 
0.50 


0.25 
0.46 


0.74 
0.52 


0.26 
0.48 


9.13 

3.88 


9.13 
3.92 


0.86 
0.91 



TABLE II: Dissipation Dk, residence probability Pr^, on k — {1,2} corresponding to A'^i LHl and N2 LII2 complexes respec- 
tively. Stoichiometry s and efficiency rj are also shown. 



Hence, for the present discussion, the most important finding from these numerical simulations is that the adaptation 
of purple bacteria does not lie in the single excitation kinetics. In particular, LLI membranes are seen to reduce their 
efficiency globally at the point where photons are becoming scarcer - hence the answer to adaptation must lie in some 
more fundamental trade-off (as we will later show explicitly). Due to the dissimilar timescales between millisecond 
absorption from eq.S and nanosecond dissipation, multiple excitation dynamics are also unlikely to occur within a 
membrane. However, multiple excitation dynamics cannot be regarded a priori not to be a reason for purple bacteria 
adaptation. A numerical study in Ref.[S] involving blockade in which two excitations can not occupy the same site and 
annihilation, where two excitations annihilate due to vibrational recombination occurring due to significant Frank- 
Condon overlap on higher exciton states, shows that these mechanisms decrease the efficiency of membranes equally 
on LLI and HLI membranes, therefore, keeping the best performance to HLI. 

Nevertheless a more comprehensive study to be presented in the next section taking into account the coupled 
exciton-RC cycling is able to justify nature's choice of proceeding with the observed bacterial adaptations. 
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FIG. 13: (a) Monte Carlo calculation of efficiency rj of HLI (diamonds) and LLI (crosses) grown membranes, as a function of 
the RC-cycling time r. Continuous lines give the result of the analytical model, (b), (c) and (d) show the distributions p{No) 
of the number of open RCs for the times shown with arrows in the main plot for HLI (filled bars) and LLI (white bars) . 



3. Complete chromatophores: Exciton and RC cycling coupled dynamics. 



It has been here repeatedly commented that RCs perform a cycle that provides the required exciton-chemical 
energy conversion. We now explain that the answer as to how adaptation can prefer the empirically observed HLI and 
LLI structures under different illumination conditions, lies in the interplay between the excitation kinetics and RC 
cycling dynamics. By virtue of quinones-quinol and cytochrome charge carriers, the RC dynamics features a 'dead' 
(or equivalently 'busy') time interval during which quinol is produced, removed and then a new quinone becomes 
available [IHIIH]. Once quinol is produced, it leaves the RC and a new quinone becomes attached. These dynamics 
are introduced into the simulation algorithm as presented in section R by closing an RC for a (random) poissonian 
distributed time with mean r after two excitations form quinol. The cycle can be sketched as follows: open RC^- 
2 ionizing excitations form quinol -^ closed RC in a time with mean r— > open RC. This RC cycling time r implies 
that at any given time, not all RCs are available for turning the electronic excitation into a useful charge separation. 
Therefore, the number of useful RCs decreases with increasing r. Too many excitations will rapidly close RCs, 
implying that any subsequently available nearby excitation will tend to wander along the membrane and eventually 
be dissipated - hence reducing rj. For the configurations resembling the empirical architectures (Fig 12), this effect 
is shown as a function of r in Fig. 13 [a) yielding a wide range of RC-cycling times at which LLI membrane is 
more efhcient than HLI. Interestingly, this range corresponds to the measured time-scale for r of milliseconds [26], 
and supports the suggestion that bacteria improve their performance in LLI conditions by enhancing quinone-quinol 
charge carrier dynamics as opposed to manipulating exciton transfer. As mentioned, a recent proposal [35, has shown 
numerically that the formation of LH2 para-crystalline domains produces a clustering trend of LHl complexes with 
enhanced quinone availability - a fact that would reduce the RC cycling time. 

However, the crossover of efficiency at r « 3 ms implies that even if no enhanced RC-cycling occurs, the HLI will 
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be less efBcient than the LLI membranes on the observed r time-scale. The explanation is quantitatively related to 



the number No of open RCs. Figs. 13 |b), (c) and (d) present the distribution p{No) of open RCs, for both HLI and 



LLI membranes and for the times shown with arrows in Fig 13 ^a). When the RC-cycling is of no importance (Fig. 
|13[b)) almost all RCs remain open, thereby making the HLI membrane more efficient than LLI since having more 



(open) RCs induces a higher probability for special pair oxidation. Near the crossover in Fig. 13 both membranes 
have distributions p{No) centered around the same value (Fig. [Tsj^c)), indicating that although more RCs are present 
in HLI vesicles, they are more frequently closed due to the ten fold light intensity difference, as compared to LLI 
conditions. Higher values of r (Fig. [Tsj^d)) present distributions where the LLI has more open RCs, in order to yield 
a better performance when photons are scarcer. Note that distributions become wider when RC cycling is increased, 
reflecting the mean-variance correspondence of Poissonian statistics used for simulation of r. Therefore the trade-off 
between RC-cycling, the actual number of RCs and the light intensity, determines the number of open RCs and hence 
the performance of a given photosynthetic vesicle architecture (i.e. HLI versus LLI). 

Hence, even though these adaptations show such distinct features in the experimentally relevant regimes for the 
RC-cycling time and illumination intensity magnitude 26l[2Hll36], FigsIIsFc) and (d) show that the distributions of 
open RCs actually overlap implying that despite such differences in growing environments, due to the adaptations 
arising, the resulting dynamics of the membranes become quite similar. Growth conditions generate adaptations that 
allow on LLI membrane to have a larger number of open RCs than the HLI adaptation and therefore LLI membrane 
will perform better than HLI with respect to RC ionization irrespective of any funneling dynamics. The inclusion of 
RC dynamics implies that the absorbed excitation will not find all RCs available, and somehow funneling would limit 
the chance of a necessary membrane sampling to explore further open RCs. Globally, a given amount of closed RCs 
will eventually alter the excitation's fate since probable states of oxidization are readily reduced. In a given lifetime, 
an excitation will find (depending on r and current light intensity I) a number of available RCs - the effective 
stoichiometry - which is different from the actual number reported by Atomic Force Microscopy [31 1361 . 

The effect of incident light intensity variations relative to the light intensity during growth with both membranes, 
presents a similar behavior. In Ref. O [6] such study is performed and it is concluded that: the LLI membrane 
performance starts to diminish well beyond the growth light intensity, while the HLI adaptation starts diminishing 
just above its growth intensity due to rapid RC closing that induce increased dissipation. Hence, in LLI membranes 
excess photons are readily used for bacterial metabolism, and HLI membranes exploit dissipation in order to limit the 
number of processed excitations. In the same work, it is found that the effect of the arrangement itself is lost due to 
RC dynamics, since the effective stoichiometry with spread out open RCs becomes alike among different membranes 
sharing the same amount of RCs, equal cycling time t and incident light intensity. 

To summarize so far, the arrangement of complexes changes slightly the efficiency of the membranes when no RC 
dynamics is included - but with RC dynamics, the most important feature is the number of open RCs. Although 
the longer RC closing times make membranes more prone to dissipation and decreased efficiency, it also makes the 
architecture less relevant for the overall dynamics. The relevant network architecture instead becomes the dynamical 
one including sparse open RCs, not the static geometrical one involving the actual clustered RCs. 

4. Analytical model 

Within a typical fluorescence lifetime of 1 ns, a single excitation has travelled hundreds of sites and explored the 
available RCs globally. The actual arrangement or architecture of the complexes seems not to influence the excitation's 
fate, since the light intensity and RC cycling determine the number of open RCs and the availability for P oxidation. 

4.1 Excitation transfer- RC cycling rate model 

Here, we present an alternative rate model which is inspired by the findings of the numerical simulations, but which 
(1) globally describes the excitation dynamics and RC cycling, (2) leads to analytical expressions for the efficiency of 
the membrane and the rate of quinol production, and (3) sheds light on the trade-off between RC-cycling and exciton 
dynamics [51 [5] . 

Shortly, Ne excitations area absorbed by the membrane at a rate 7^1, and will find its way to become RC ionizations 
with a rate per particle Xc{No) whose dependence on the number of open RCs is made explicit, or be dissipated at a 
rate 713. On the other hand, RCs have their own dynamycs, closing at a rate Xc{No) x Ne and individually opening 



at a rate 1/t. The dependance Xc{No) is numerically found aidded by the stochastic model and shown in Fig 14 ^a 
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FIG. 14: (a) Numerical results showing the rate of ionization Xc{No) of an RC for HLI (diamonds) and LLI (crosses) membranes, 
together with a quadratic (dashed line) and linear (continuous) dependence on the number of closed RCs (TVi — No). The fitting 
parameters for a + bNo are a = {15.16, 7.72}ns"\ b = {-0.21,-0.21} ns"^ and for a + bNo + cN^, a = {16.61, 8.21}ns"\ 
b — {—0.35, —0.33}, and c — {3.6, 1.5}/is~^, for HLI and LLI membranes respectively, (b) rj as function of r and a = I/Io, 
obtained from the complete analytical solution for LLI (white) and HLI (grey) membranes 



Rate equations can therefore be written: 



dNE 

dt 

dNo 

dt 



= -iXc{No)+lD)NE+lA 

r / 



to be solved, and of use for the calculation of the steady-state efhciency 77 — njic/nA- 

Xc{No)Ne 
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lA 



(23) 
(24) 

(25) 

(26) 

(27) 



where A^ is the rate of RC ionization when no RC-cycling is accounted, dependent only on the amount of RCs 



A linear fit for \c(No) allows an analytical expression for r\ 

1 



r7(r, 7^(/)) = -^-^ {2iVi(A°, + 713) + 7aA^i 
^^Nl{\% + ^n? + 4iVi7AA^(7D - \1)t + {lAXlrf 



present in the vesicle jS]. This analytical expression is shown in Fig|l4[b) and illustrates that 77 > 0.9 if the transfer-P 

As can be seen in Figs. 
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reduction time is less than a tenth of the dissipation time, not including RC cycling, 
the analytical solution is in good quantitative agreement with the numerical stochastic simulation, and provides 
support for the assumptions made. Moreover, this model shows directly that the efficiency is driven by the interplay 
between the RC cycling time and light intensity. Figure |14[ b) shows up an entire region of parameter space where 
LLI membranes are better than HLI in terms of causing P ionization, even though the actual number of RCs that 
they have is smaller. In view of these results, it is interesting to note how clever Nature has been in tinkering with 
the efficiency of LLI vesicles and the dissipative behavior of HLI adaptation, in order to meet the needs of bacteria 
subject to the illumination conditions of the growing environment. 



4.2 Bacterial metabolic demands 



Photosynthetic membranes must provide enough energy to fulfill the metabolic requirements of the living bacteria 
quantified by the quinol output or quinol rate 
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FIG. 15: (a) Quinol rate W in HLI (diamonds, Iq — lOOW/ni^) and LLI (crosses, Iq — lOW/m^) grown membranes, as a 
function of incident intensity I. (b) Quinol rate contours oi W — {1900,2000,2100,2200} s^^ in black, blue, red and pink, 
respectively. 



which depends directly on the excitations that ionize RCs fine- The factor ^ accounts for the requirement of two 
ionizations to form a single quinol molecule. Although these membranes were grown under continuous illumination, 
the adaptations themselves are a product of millions of years of evolution. Using RC cycling times that preserve 



quinol rate in both adaptations, different behaviors emerge when the illumination intensity is varied (see Fig. 15 'a). 
The increased illumination is readily used by the LLI adaptation, in order to profit from excess excitations in an 
otherwise low productivity regime. On the other hand, the HLI membrane maintains the quinone rate constant, 
thereby avoiding the risk of pH imbalance in the event that the light intensity suddenly increased. We stress that 
the number of RCs synthesized does not directly reflect the number of available states of ionization in the membrane. 
LLI synthesizes a small amount of RCs in order to enhance quinone diffusion, such that excess light intensity is 
utilized by the majority of special pairs. In HLI, the synthesis of more LHl-RC complexes slows down RC-cycling, 
which ensures that many of these RCs are unavailable and hence be advantageous of evenly distributed dissipation to 
steadily supply quinol independent of any excitation increase. The very good agreement between our analytic results 
and the stochastic simulations, yields additional physical insight concerning the stoichiometrics found experimentally 
in Rsp. Photometricum. 

A closed form expression regarding all dynamical parameters involved can be obtained [S]: 
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where the dependence on stoichiometry is made explicit due to absoprtion cross section in ja and on B{s) — 

ing fraction 

and t(s) is 



^ f(s)A 1 "^liidi is a parameter that depends on area ^1(2) of individual LHl (LH2) complexes and filling fraction 
/(s). The filling chromatophore fraction dependence on s is available from experimental data of Ref. | 
constructed from an interpolation scheme [5J. 

As emphasized in Ref. [34] , membranes with s~Q or s~2 were not observed, which is to be compared with the 
contour plots regarding constant quinol output, of s as a function of growing light intensity /q, shown in Fig |15[ b). 
These results support a dichotomic observation where s « 4 predominantly on a great range for growing light intensity. 
However, a prediction can be made for 30-40 W/m^ where a great sensitivity of stoichiometry ratios rapidly build up 
the number of antenna LH2 complexes. Very recently [25j . membranes were grown with 30W/m^ and an experimental 
stoichiometry of 4.8 was found. The contour of 2200 s~^ predicts a value for stoichiometry of 4.72 at such growing light 
intensities. This agreement is quite remarkable, since a simple linear interpolation among the values s w 4, /q = 100 
W/m^ and s w 8, /q = 10 W/m^ would wrongly predict s = 7.1 at 30 W/m^. We encourage experimentalists to 
confirm the full range of predicted behaviors as a function of light-intensity and stoichiometry. Such exercise would 
without doubt confirm/expand the understanding on RC-exciton dynamics trade off, pinpointing a direction to pursue 
solar energy conversion research, provided by Nature itself. 
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4.3 Performance of photosynthetic membranes under extreme photon statistics 

Photosynthetic (e.g purple) bacteria provide the crucial coupling between the Sun's energy and the production of 
food on Earth, and have adapted successfully to a variety of terrestrial conditions since the beginnings of life on Earth 
several billion years ago. In this section we explore whether terrestrial bacteria, which are the product of millions of 
years of evolutionary pressure on Earth, could survive if suddenly exposed to incident light with extreme statistics. 
We are therefore mimicking a scenario in which purple bacteria are either (i) suddenly transported to some unknown 
extreme solar environment elsewhere in the universe, or (ii) where our own Sun suddenly picks up extremal behavior 
in terms of the temporal statistics of its emitted photons, or (iii) the bacteria are subjected to extreme artificial 
light sources such as that in Ref. [4] involving stimulated Raman scattering, as well as in coherent anti-Stokes Raman 
scattering in silicon, or (iv) the external membrane skin of the bacterium is modified in such a way that the absorption 
of photons takes on extreme statistical properties. 

Our discussion is qualitatively different from previous work looking at life in extreme conditions, since those dis- 
cussions have tended to focus on environmental extremes affecting the biochemical metabolism or extremes in the 
incident light intensity. By contrast, our current discussion focuses on extremes in the incident photon statistics. 
Although the setting is largely hypothetical, our theoretical predictions are based on a realistic semi-empirical model 
which incorporates (i) high precision empirical AFM information about spatial locations of LHl and LH2 biomolec- 
ular complexes in the membrane architecture of Rsp. Photometricum purple bacteria, and (ii) full-scale stochastic 
simulations of the excitation kinetics and reaction center dynamics within the empirical membrane. 

As with any process involving events occurring in a stochastic way over time, the statistical properties of arriving 
photons may show deviations from a pure coin-toss process in two broad ways: burstiness and memory jl5| . First 
consider the simplest process in which the rate of arrival of a photon has a constant probability per unit time. It 
is well known that this so-called Poisson process produces a distribution of the waiting time for the next photon 
arrival which is exponential in form, given by Pp{t) ~ exp(— r/rg). The extent to which the observed arrival time 
distribution P(t) deviates from exponential, indicates how non-Poissonian the photon arrival is. Following Barabasi 
[Is] , we refer to this as 'burstiness' B and define it by its deviation from a purely Poisson process: 

^ {(Tr/mr - 1) _ (o-^ - TUr) 
[Gr/mr + 1) {(Tt + nir) 

where CTt and tti,- are the standard deviation and mean respectively of the empirical distribution P{t). For a pure 
Poisson process, the mean and standard deviation of the arrival time distribution are equal and hence B — 0. The 
other property which can be noticeable for a non-Poisson process is the memory M between consecutive inter-arrival 
times which, following Barabasi, we define as: 

j^j^ " ^^ (t, -mi)(T,+i -TO2) 






Ur - i ^ 0-lfT2 

where Ut is the number of inter-arrival times measured from the signal and mi(m2) and <ti(i72) are sample mean and 
sample standard deviation of Tj(ri+i)'s respectively {i ~ 1, . . . ,nr — !)■ For a pure Poisson process, M = 0. Both 
B and M range from —1 to +1, with the value B — and M = for a strict Poisson process. We will assume for 
simplicity that the arriving photons are all absorbed and hence one exciton created within the membrane by each 
photon. This can easily be generalized but at the expense of adding another layer of statistical analysis to connect 
the statistics of arriving photons to the statistics of the excitons being created within the LII2/LII1 membrane - 
indeed, taking a constant absorption probability less than unity would not change our main conclusions. We will also 
neglect the possibility that several photons arrive at exactly the same time. In principle, incident photons can be 
generated numerically with values of B and M which are arbitrarily close to any specific (M, B) value - however this 
is extremely time-consuming numerically. Instead, we focus here on specific processes where the B and M can be 
calculated analytically. Although this means that the entire (M, B) parameter space is not accessed, most of it can 
indeed be - and with the added advantage that analytic values for B and M are generated. 

Figure U& summarizes our findings in terms of the incident photons (i.e. excitation input) and metabolic output 
from the LH2/LII1 membrane, over the entire (M, B) parameter space. The subspace shown by the combination of 
the white and dotted regions (i.e. the region in the (M, B) space which is not diagonally shaded) comprises points 
with (Af, B) values that can be generated by one of three relatively simple types of photon input: (a) step input, (b) 
bunched input and (c) power-law step input, as shown in the three panels respectively. Each point in the blank or 
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dotted region denotes a time-series of initial excitations in the membrane with those particular burstiness and memory 
values (-B and M). This train of initial excitations then migrates within the membrane of LH2/LH1 complexes, subject 
to the dynamical interplay of migration and trapping as discussed earlier in this paper, and gives rise to a given output 



time-series of 'food' (quinol) to the bacteria. Figure 16 shows explicitly a variety of initial (M, B) values (crosses 
and stars) and the trajectory represents the locus of resulting quinol outputs from the reaction centers (RC) for this 
particular (M, B) input. The trajectory is generated by varying the RC closing time within the range of physically 
reasonable values. The trajectory is finite since the range of physically reasonable RC closing times is also finite (20 
to lOOOHz). As the time during which the RC is closed increases, the output becomes more Poisson-like, i.e. the B 
and M output values from the membrane after absorption at the RC, tend towards 0. Hence the trajectories head 
toward the center as shown. 

The (M, B) values for natural sunlight on Earth lie near M = and B = since the incident photon arrival from 
the Sun is approximately a Poisson process. This produces a quinol output which is also a Poisson process (M = 
and -8 = 0). Now suppose a terrestrial bacteria is suddenly subjected to an extreme incident light source with (M, B) 
values which lie at a general point in (M, B) parameter space. Can it survive? It is reasonable to expect that the 
bacteria will not survive if the resulting quinol output is very different from that on Earth (i.e. very different from 
B — and M = 0) since the bacteria is well-tuned for Life on Earth only. Hence we say that the bacteria can only 
survive if the quinol output has (M, B) values within 0.01 of the terrestrial values of B = and M = 0. The blank 
spaces denote (M, B) values for which the bacteria would survive, while the dotted spaces are where the bacteria 
would die. 

Remarkably, there are therefore substantial regions of extreme incident photon statistics (i.e. B y^ and M ^ 0) 
where the bacteria would survive, producing a quinol output for the rest of the organism which resembles that on 
Earth (i.e. \B\ < e and \M\ < e where e = 0.01). The fact that these blank spaces have an irregular form is due 
to the highly nonlinear nature of the kinetic interplay between exciton migration and trapping, as discussed earlier 
in this chapter. For example, panel (c) shows that a bacterium which is only tuned to survive on Earth can have a 
maximally negative burstiness B = —1 and a maximal memory M — 1, meaning that the incident photons arrive like 
the regular ticks of a clock. The same applies to all the incident photon conditions demonstrated by the blank spaces. 
If we make the survival criterion more generous (i.e. e > 0.01) the blank spaces corresponding to survivability of the 



terrestrial bacterium increase in size until no dotted space remains (i.e. when e = 1). All the dotted areas in Fig. 16 
become blank spaces and hence represent survivable incident photon conditions for the bacteria. 

5. Perspectives: Photosynthetic membranes of purple bacteria as a basis to develop stable energy conversion 

This chapter has covered a comprehensive review of purple bacteria adaptation to light intensity conditions and has 
provided a basis to understand the chromatic adaptation of photosynthetic harvesting membranes as a consequence of 
the requirement imposed by nature on living organisms to develop a machinery capable of fuel production in a stable 
manner. It is worth to emulate, and certainly interesting to develop in artificial devices to surpass the usual bleaching 
on semiconductor panels, an analogous mechanism to the excitation kinetics and RC dynamics interplay, in order to 
shift the emphasis of requirements of fuel production and allow highly efficient energy transfer to available charge 
carriers (quinones) in low light intensity conditions, while consistently, a leveled fuel (quinol) production statistics is 
accomplished when photon arrival varies both in intensity and photon waiting time statistics. 
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